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ABSTRACT 

We present CCD photometry in UBVRI of several thousand Galactic field stars in 
four large (> 1 degree 2 ) regions centered on diffuse interstellar dust clouds, commonly 
referred to as "cirrus" clouds (with optical depth Ay less than unity). Our goal in 
studying these stars is to investigate the properties of the cirrus clouds. A comparison of 
the observed stellar surface density between on-cloud and off-cloud regions as a function 
of apparent magnitude in each of the five bands effectively yields a measure of the 
extinction through each cloud. For two of the cirrus clouds, this method is used to 
derive UBVRI star counts-based extinction curves, and [/-band counts are used to place 
constraints on the cloud distance. The color distribution of stars and their location in 
(U — B, B — V) and (B — V, V — I) color-color space are analyzed in order to determine 
the amount of selective extinction (reddening) caused by the cirrus. The color excesses, 
A\ — Ay, derived from stellar color histogram offsets for the four clouds, are better fit 
by a reddening law that rises steeply towards short wavelengths [Ry = A V /E(B — V) ^ 
2] than by the standard law (Ry = 3.1). This may be indicative of a higher-than- 
average abundance of small dust grains relative to larger grains in diffuse cirrus clouds. 
The shape of the counts-based effective extinction curve and a comparison of different 
estimates of the dust optical depth (extinction optical depth derived from background 
star counts/colors; emission optical depth derived from far infrared measurements), are 
used to measure the degree of clumpiness in clouds. The set of techniques explored in 
this paper can be readily adapted to the Sloan Digital Sky Survey data set in order to 
carry out a systematic, large-scale study of cirrus clouds. 

Subject headings: ISM: clouds — dust, extinction — ISM: general — ISM: structure — stars: 
fundamental parameters — methods: statistical 
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1. Introduction 

Among the most remarkable of the many discov- 
eries made by the Infrared Astronomical Satellite 
(IRAS, Neugebauer et al. 1984) was the detection of 
a bright component of the far infrared background 
emission. This far infrared component is marked 
by discrete condensations of filamentary clouds with 
characteristic sizes ranging between a few and ~ 30°. 
Their resemblance to clouds in our own atmosphere 
led Low et al. (1984) to dub them the infrared "cir- 
rus". The cirrus clouds are most readily observed at 
high Galactic latitudes where there is relatively little 
contribution to the net background by the Galactic 
plane, and are most prominent at 100 /im, although 
they are visible in all four IRAS bands. Low et al. 
noted the positional correlation of some of the cir- 
rus emission with Hi clouds mapped by Heiles (1975) 
and interpreted this as evidence that the infrared cir- 
rus was in fact associated with the diffuse interstellar 
medium (ISM). 

Despite the fact that optical cirrus was first dis- 
covered over four decades ago (de Vaucouleurs 1955, 
1960), definite association with high latitude dust 
came much later (de Vaucouleurs & Freeman 1972; 
Sandage 1976). The first systematic comparison be- 
tween the IRAS 100 /im and optical morphologies of 
dust clouds was conducted by de Vries & Le Poole 
(1985). This and a few other studies (Paley 1990; 
Stark 1993), all based on photographic data, reveal 
a good correspondence between optical and infrared 
cirrus. Quantitative optical surface brightness mea- 
surements of cirrus clouds have been carried out by 
Guhathakurta & Tyson (1989), Paley et al. (1991), 
and Gordon, Witt, & Friedmann (1998), while spec- 
troscopy of diffuse cirrus clouds has been carried 
out only recently (Szomoru & Guhathakurta 1998). 
Such optical measurements provide important clues 
towards understanding grain optical properties, such 
as the albedo and the phase function asymmetry, 
and have demonstrated the existence of extended red 
emission, caused by photoluminscence in very small 
grains, in the diffuse ISM (Szomoru & Guhathakurta 
1998; Gordon et al. 1998). 

Diffuse cirrus clouds typically extend over large ar- 
eas on the sky. Such clouds can be probed through 
the study of their effect on the light of background 
sources (Galactic stars, distant quasars). A statistical 
method which is well suited for cirrus clouds employs 
star counts to determine the extinction. This method, 



which was first described and applied several decades 
ago (Wolf 1923; Bok 1937), has been used to deter- 
mine the extinction in high latitude molecular clouds 
(cf. Magnani & de Vries 1986; Stark 1995). By com- 
paring cumulative star counts in the direction of an 
interstellar cloud with star counts in extinction-free 
regions, or alternatively with model-based star counts 
(so-called Wolf diagrams — Wolf 1923), one may de- 
rive both the amount of extinction through and the 
distance to the cloud. An alternate method relies 
on optical and ultraviolet spectroscopy of individual 
background hot stars of known spectral type (Stark 
1995) or QSOs with mostly featureless, power-law 
spectra (cf. Bowen 1991; Bowen, Blades, & Pettini 
1995). This latter method is useful for the study of 
dust extinction curves, cloud chemical abundances, 
and excitation conditions in the gaseous ISM. It has 
been applied mostly to dense molecular clouds cores 
with Ay y> 1 and to the warm ionized ISM of the 
Galaxy, rather than to cirrus clouds. 

This paper presents photometry in UBVRI of sev- 
eral thousand stars in four large area (> 1 degree 2 ) 
fields centered on cirrus clouds. The surface photom- 
etry of these clouds will be the subject of a future pa- 
per. Our study differs from previous studies in several 
respects: 

• CCD photometry of faint stars over a wide field 
is less subject to systematic error than photo- 
graphic measurements thanks in large part to 
accurate flat fielding and sky subtraction. 

• The use of UBVRI bands opens up the possibil- 
ity of studying the effect of selective extinction 
by interstellar dust on stellar color distributions 
and color-color diagrams, in addition to its ef- 
fect on star counts. This combination of tech- 
niques can be used to investigate the occurence 
of dense cores and to quantify t he amount of 
structure within the clouds (Sec. [T^). 

• Measurements of the effect of extinction on star 
counts/color distribution are made relative to 
a well-matched off-cirrus sample of stars, and 
this yields a reliable and precise determination 
of the extinction curve (A\). 

• Data at short wavelengths (3600A U band) en- 
ables us to study clouds of relatively low optical 
depth; all four clouds studied here have Ay val- 
ues significantly less than unity. Most previous 
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studies employing background star counts have 
been of molecular clouds with higher optical ex- 
tinction than the cirrus clouds in our sample, 
with the study by Stark (1995) being a notable 
exception; in fact, two of our four clouds are also 
in his sample. In contrast to our UBVRI CCD 
data, however, the optical part of his study is 
based solely on £/-band photographic plates. 

Direct measurements of the optical extinction have 
a number of important applications. Firstly, the ratio 
of Ay (visual extinction optical depth) to rioo^m (the 
far infrared emission optical depth) provides a test, 
albeit a crude one, of interstellar dust grain mod- 
els. Such tests are now possible thanks to accurate 
Ay measurements, along with recent improvements 
in techniques for the processing of far infrared data 
and hence in the measurement of rioo^m (Schlegel, 
Finkbeiner, & Davis 1998, hereafter SFD; Verter & 
Rickard 1998; Verter et al. 1998). Secondly, a di- 
rect measurement of Ay (as opposed to an approxi- 
mate scaling based on the far infrared surface bright- 
ness) is essential in order to model the degree of pen- 
etration and internal reddening of ambient starlight 
in cirrus clouds; this is a key piece of information 
used in the interpretation of ultraviolet /optical/near- 
infrared spectra and U BVRI surface brightness mea- 
surements of optical cirrus. 

The observations and reduction of optical data are 
described in Sec. Hand far infrared measurements are 
described in Sec. [| Sec. || contains the background 
star count analysis along with the resulting determi- 
nation of the extinction law and cloud distance; Sec. ^ 
contains a description of stellar color-color diagrams 
and their use in the measurement of selective extinc- 
tion (reddening). In Sec. ^, the optical extinction is 
compared to the 100 /im optical depth, and optical 
depth estimates from the various methods are dis- 
cussed in the context of cloud structure (i.e., dumpi- 
ness) and dust properties. Sec. [7] contains a summary 
of the main results. 

2. Optical Data 

2.1. Targets and Observations 

Over the last decade, several compact, relatively 
isolated, high latitude cirrus clouds have been ob- 
served at optical wavelengths for the purpose of study- 
ing the reprocessing of starlight by interstellar dust 
grains (cf. Guhathakurta & Cutri 1994). These cirrus 



targets have been selected from IRAS 100 fim maps 
or drawn from optical cirrus catalogs such as the list 
of Lynds Bright Nebulae (Lynds 1965) or the compi- 
lation of Paley (1990) which are, in turn, derived from 
Palomar Observatory Sky Survey plates. This paper 
is based on observations of four such cirrus clouds, 
RCrA (R Corona Australis), PV1 (Project Verifica- 
tion #1, an early-phase target of the Infrared Space 
Observatory), Paley 1, and Paley 3 (from the Paley 
1990 catalog). 

The observations were carried out with the Cerro 
Tololo Interamerican Observatory^ Curtis-Schmidt tele- 
scope, equipped with a Tektronix 2048 x 2048 CCD 
camera, during three nights in 1995 October. This 
telescope has a 0.9-m primary mirror and a 0.6-m 
corrector. The pixel size of 21 /im, corresponding 
to a scale of 2 '.'03 pixel -1 , yields a field of view of 
1.15° x 1.15°. The finite filter size causes slight vi- 
gnetting (~ 20% at the extreme corners of each CCD 
frame) but this is well corrected by sky flat fields. 

The Schmidt CCD observations consist of a se- 
ries of disregistered 400 s exposures in the standard 
UBVRI filter set at Cerro Tololo Interamerican Ob- 
servatory. The effective wavelengths of these filters 
are practically identical to those of the corresponding 
Johnson bandpasses: 3597A, 4405A, 5495A, 6993A, 
and 9009A, respectively. The total integration time 
for RCrA, PV1, and Paley 3 is about 20-30 min per 
field in each of the BVRI bands, and about 50 min in 
the U band where the instrumental efficiency is very 
low. The integration times for the Paley 1 field are 
about a factor of 2 longer than in the corresponding 
bands for other three fields. For photometric calibra- 
tion purposes, the Landolt Selected Area SA 92 (Lan- 
dolt 1992) containing several standard stars was ob- 
served in UBVRI twice during the first night. Short 
twilight flat exposures in all five bands were obtained 
at the beginning and end of each night; several se- 
ries of bias frames were obtained during the daytime. 
The FWHM of the stellar images ranges from 5'.'1- 
5 '.'7 (2.5-2.8 pixels), the result of atmospheric seeing, 
imperfect focus, and coarse pixel scale. The coordi- 
nates and observational parameters of the clouds are 
listed in Table [jj 



3 Cerro Tololo Interamerican Observatory, National Optical As- 
tronomy Observatories, is operated by the Association of Uni- 
versities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation. 
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Table 1 

Cloud Coordinates and Observational Parameters 
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tint (U) 
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tint(I) 




m 


(° ' ") 
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(°) 


(min) 


(min) (min) (min) 


(min) 



RCrA 


18 56 04.1 


-37 18 05 


359.2 


-16.9 


47 


27 


20 


20 


27 


PV1 


22 06 28.4 


-03 32 33 


56.5 


-44.1 


53 


27 


20 


20 


27 


Paley 1 


02 38 24.6 


-29 43 36 


225.6 


-66.4 


93 


53 


47 


40 


47 


Paley3 


23 48 53.8 


-71 47 12 


309.7 


-44.5 


53 


33 


20 


20 


33 



Note. — 

Col. 6 - 10: Total integration time in minutes in each of the U, B, V, R, and I bands. 



2.2. Data Reduction 

As a first step, all the images (object, photomet- 
ric calibration, and flat-field exposures) are overscan- 
and bias-subtracted and trimmed. The twilight flats 
and selected disregistered (dark sky) object exposures 
of all four cirrus clouds are combined into a master 
flat-field image for each of the UBVRI bands; the 
master fiat-field image is applied to all object and 
photometric calibration exposures. After flat field- 
ing, image defects (e.g., hot pixels, charge traps, and 
bad columns) are replaced by the median value of the 
pixels surrounding the defects. 

A set of relatively bright and isolated stars are 
selected in each of the four cirrus fields, and these 
are used to determine the linear shift in (x, y) be- 
tween each of the images and a particular astromet- 
ric reference image. After registering the images ap- 
proximately, the IRAF routine GEOMAP is applied to 
these same bright stars in order to determine higher 
order coordinate transformations (translation, rota- 
tion, magnification, distortion) for every image with 
respect to the reference image; the IRAF routine GEO- 
tran is then used to accurately align all the images 
using the coordinate transformations determined by 

GEOMAP. 

The median background sky level for each object 
exposure is subtracted from the image. The regis- 
tered, sky-subtracted images of a given cirrus cloud 



in a given band are then combined into a median- 
averaged image. In the case of the PV 1 field, the 
translational offsets between various exposures are 
quite large (of the order of the field of view of an 
individual CCD image); the final combined image in 
each band for PV 1 is the union of all the images in 
that band. As a result, the effective exposure time, 
and hence the r.m.s. sky noise level, varies with posi- 
tion across the final image of the PV 1 field. For the 
other three cirrus fields, only the intersection area of 
all exposures is used, and the noisy (partial overlap) 
edges of the combined images are discarded. 

The final, combined V-band images of the four cir- 
rus cloud fields are shown in Figure [l]. The images 
contain several very bright, saturated stars with ex- 
tensive and complicated charge bleed patterns; the 
exact number of saturated stars in the image depends 
on the Galactic latitude of the field (most numerous 
in RCrA) and on wavelength (many more in R- and 
J-band images than in UBV). Instead of attempting 
to model (and subtract) the intensity distribution of 
saturated stars, the affected regions in their vicinity 
are eliminated from further consideration, since de- 
tection and photometry of faint stars is unreliable in 
these regions. This is done by finding the positions of 
stars for which the central pixel value is above the sat- 
uration limit (~ 60, 000 ADU or 1.26 x 10 5 electrons) 
and determining the number of saturated pixels at 
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the center of every saturated star. An empirical rela- 
tion is then defined between the number of saturated 
pixels at the center of a stellar image and the size of 
the associated region to be eliminated; this relation 
is used to mask a circular area of appropriate radius 
around every saturated star. 

2.3. Stellar Photometry 

Lists of object positions are generated indepen- 
dently for each of the five bands and for each of the 
four cirrus fields using the peak finding algorithm 
find of the stellar photometry program daophot 
(Stetson 1987, 1992). An object detection thresh- 
old of 3.5cr is used, where a is the sky noise; in the 
case of the PV 1 field, the variation of sky noise level 
across the image (due to non-uniform effective expo- 
sure time) is taken into account. 

Most of the detected objects are Galactic stars; 
a small fraction are distant field galaxies. Even at 
high Galactic latitudes > 30°), the surface den- 
sity of stars exceeds that of field galaxies for V <J 20 
(cf. Kron 1980), roughly the apparent magnitude limit 
of our study (see Sec. 2.4 and Figs. || and ||). Given 
the relatively coarse angular resolution of our CCD 
data and the fact that the majority of the galaxies 
in the sample are expected to be close to the detec- 
tion/completeness limit, no attempt has been made 
to exclude galaxies using morphological star-galaxy 
separation. In fact, field galaxies are just as useful 
as distant stars for studying the optical depth of a 
foreground cirrus cloud. In the rest of this paper, we 
loosely use the term "stars" to refer to the sample of 
detected (mostly stellar) objects. 

An iterative procedure, based on a combination of 
several daophot tasks, is applied to the daophot/find 
star list in order to build an empirical point spread 
function template and to fit the template to every 
star in the list. Aperture photometry is carried out 
for each star in turn, with the results of the point 
spread function fit used to subtract off the light of 
the neighbors of the star in question. This combi- 
nation of point spread function fitting and aperture 
photometry techniques is ideal for photometry of un- 
dersampled stellar images in crowded fields (Guha- 
thakurta et al. 1996). The aperture magnitudes, cal- 
culated over small (2 pixel radius) apertures to avoid 
neighbor contamination, are corrected to total mag- 
nitudes using standard "curve-of-growth" corrections. 
The curve of growth is determined independently for 
each band and for each cirrus field, based on the in- 



tensity profiles of a set of isolated bright stars. The 
lists of star positions and corresponding total instru- 
mental magnitudes are then combined into a position- 
matched UBVRI list for each cirrus field. Total in- 
strumental magnitudes are also determined for pho- 
tometric standard stars in the SA 92 calibration field 
using a similar curve-of-growth method. 

As the SA 92 calibration field was only observed 
on the first (photometric) night, a bootstrap mag- 
nitude zeropoint adjustment is determined for each 
field and for each band: the airmass-corrected instru- 
mental magnitudes of a sample of secondary stan- 
dard stars in a single image obtained during that 
first night are compared to those measured on the 
median-averaged image. This zeropoint adjustment 
is quite small, ranging from 0.01 - 0.09 mag. One of 
the four fields (Paley3) was not observed during the 
first night and thus could not be bootstrap corrected. 
However, the resulting error in the overall magnitude 
scale for this field is likely to be small (< 0.1 mag) 
and thus unimportant for the analysis in this paper. 

The IRAF routine fitparam is used to derive the 
zeropoints and color- and airmass-coefficients of the 
transformation equations for converting instrumen- 
tal magnitudes to calibrated Johnson-Kron-Cousins 
UBVRI magnitudes, based on measurements of the 
total instrumental magnitudes of photometric stan- 
dard stars in the SA92 field. The instrumental stel- 
lar magnitudes in the four cirrus fields are converted 
to calibrated UBVRI magnitudes using the routine 
invertfit, adopting an airmass term that is the av- 
erage of all the individual exposures for that band 
in each field. The transformation equations from 
(UBVRI) ms t — > (UBVRI)caiib contain color terms 
based on U - B, B - V, B - V, V - R, and R - 1, re- 
spectively. The corresponding color term coefficients 
are -0.034, +0.165, -0.073, +0.027, and -0.013. Not 
all stars are detected in all bands; the combination 
of B VI detections is of particular interest to us (see 
Sec. 5.2), including a small subset of stars which hap- 
pen to be undetected in the R band. To estimate 
the color term in the L ins t — > J C aiib relation for these 
i?-band non-detections, their calibrated R — I color 
is estimated using the empirical correlation between 
(B — V^)caiib and (R — I) ca iib that is seen for stars de- 
tected in all bands. This technique is accurate enough 
for our purposes since the /-band color term coeffi- 
cient is very small. 

For the star count analysis described in Sec. |[ 
star lists are also constructed on the basis of detec- 
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tions in a single band (e.g., U band only, B band 
only, etc.). These single-band detections are cali- 
brated using an average color term in the instru- 
mental — > Johnson conversion relation; for example: 
C^caiib = CAnst + (f/calib - CAnst), where the average 
is computed over all stars for which a proper trans- 
formation is possible. Thus the calibrated UBVRI 
magnitudes in the single band lists arc only approx- 
imations to true calibrated magnitudes with proper 
color terms. This approximation however should be 
quite good as the color term coefficients are generally 
small. 

2.4. Completeness 

The completeness of the single-band star lists as 
function of magnitude, N \, s /Ntrae> is modeled as a 
double exponential: 



N ohs /N tTUe = 0.5 [2 - exp((m - m 50 )/Am)] (1) 
m<ra5o 

N ohs /N tTUe = 0.5 [exp(-(m - m 50 )/Am)] (2) 
m > mso 

The underlying star count versus apparent magnitude 
relation, A true , is assumed to be a power law in the 
vicinity of the 50% completeness point, TO50. The val- 
ues of the two free parameters 77150 and Am (sharp- 
ness with which incompleteness sets in) are then de- 
termined by fitting to the shape of the observed count 
distribution, -/V b s . At magnitudes brighter than the 
completeness limit, the observed counts are a good 
fit to the Bahcall-Soneira model (Bahcall & Soneira 
1980; see Sec. O and Fig. ||), which predicts A^rue to 
be roughly a power law over the magnitude range of 
interest. The typical value of Am is about 0.4, with 
a range from 0.3 to 0.6, while the 50% completeness 
values range from /50 w 18 to U50 ~ 21 (as indicated 
in Figs. H and |). 



in regions without visible optical cirrus. An "off- 
cirrus" region is defined in which 



F 



100/im 



< (F 



lOOjum/bg 



5(7ioo/j 



3. Far Infrared Data 

The IRAS 100 /im (ISSA) maps provided by the 
Infrared Processing and Analysis Center are used to 
divide each optical image into three regions on the ba- 
sis of 100 /im surface brightness, Fioo^m- The mean 
brightness, (Fioo^m)bg, and pixel-to-pixel r.m.s. vari- 
ation, (Xioofim, of the 100 /jm background is computed 



< (FioOum)bg + 



along with two "on-cirrus" regions: 

(FloOjum/bg + 5(7i00/jm < FloO/an ^lOOpm/bg 
HfTlOOMm ("Onl") 

and 

Fl00/im > (Fi00/im)bg + HcioO/im ("on2"). 

Each star in the matched UBVRI list and single-band 
lists is then flagged according to the region in which 
it falls. The background 100 /mi flux, r.m.s. variation, 
and the areas of the "off" , "onl" , and "on2" regions 
for each cirrus field are listed in Table |[ These areas 
are indicated on the y-band images of the four clouds 
shown in Figure [l]. 

As discussed in Sec. || below, the optical depth of 
a cirrus cloud can be estimated both from its effect 
on the counts/colors of background stars (extinction 
optical depth) and from its far infrared brightness 
and color temperature (emission optical depth). The 
far infrared emission-based estimate is derived from 
full sky 100 /im maps provided by SFD. These maps 
are a reprocessed composite of the COsmic Back- 
ground Explorer / Diffuse InfraRed Background Ex- 
periment (COBE/DIRBE) and IRAS (ISSA) maps, 
with zodiacal foreground and confirmed point sources 
removed. The procedure combines the relatively high 
angular resolution of IRAS 100 /im data (3') with 
the superior photometric calibration of DIRBE data, 
and is estimated to be twice as accurate as the older 
Burstein & Heiles (1982) reddening maps in regions 
of low to moderate reddening. In addition to 100 /im 
flux maps with improved photometric calibration, 
SFD compute Ffoou m over the full sky, the equiva- 
lent 100 /im flux corrected for variations in mean dust 
temperature from one line of sight to another. The 
amount of reddening, E(B — V), is calculated by nor- 
malizing the amplitude of reddening per unit of cor- 
rected 100 /xm flux, based on the observed redden- 
ing towards brightest cluster galaxies and elliptical 
galaxies. The E(B — V) values of SFD are scaled 
by R v = A V /E(B - V) = 3.1 to derive differential 
on-cloud A v values ("onl" -"off", "on2"-"off") for 
each of the four cirrus clouds in this study (Table || 
and Fig. H). 

SFD present all-sky maps of the dust temperature, 
averaged along the line of sight, derived from the ra- 
tio of DIRBE 100 /im to 240 /im flux measurements. 
The effective angular resolution of the temperature 
measurements is limited to ~ 1° by the resolution of 
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the DIRBE far infrared data. The mean dust tem- 
perature of each of the four cirrus clouds in this study 
(averaged over "onl" and "on2" regions) is listed in 
Table 0. The implications of cloud-to-cloud variations 



in dust temperature will be discussed in Sec. 6.3 



SFD assume a standard extinction law with Ry = 
3.1 in deriving E(B — V) values which should be cor- 
rect on average over the full sky. The lines of sight to- 
wards the four clouds in this study, however, may have 
Rv values that are different from this standard value 
(see Sec. 4.2 and 5.1). For a line of sight with a non- 
standard extinction law, the SFD value of E(B — V) 
will be inaccurate, but the estimate of Ay should be 
valid. SFD's estimate is effectively based on Tioo/jm 
(correcting the 100 /im flux to the average dust tem- 
perature). The variation in extinction law from one 
line of sight to another appears to be mostly in terms 
of the relative amount of extinction at wavelengths at 
or shortwards of the B band while the portion red- 
ward of the V band is invariant in shape (Cardelli, 
Clayton, & Mathis 1989, hereafter CCM; Martin & 
Whittet 1990). In other words, our estimate of Ay 
is valid under the assumption that extinction laws 
along different lines of sight, normalized by TiQ 0fim , 
are similar for A > Ay and diverge only at shorter 
wavelengths. 

4. Star Counts 

4.1. Analysis of Cumulative Star Counts 

Single-band star lists are used to derive star counts 
in each band and for each region ("off", "onl", and 
"on2") of each field. The single-band star lists are 
preferable to the matched U BVRI lists for this pur- 
pose, as they are based on clean, well-defined selec- 
tion functions and contain more stars by virtue of a 
higher degree of completeness at the faint end. Fig- 
ure U shows cumulative star counts versus apparent 
magnitude, also called Wolf diagrams (Wolf 1923), for 
the four cirrus fields in the U band; Figure || shows 
Wolf diagrams for RCrA in BVRI. As these dia- 
grams are based on single-band detections, the pho- 
tometric conversion to the Johnson magnitude system 
is only approximate (Sec. ||). The "onl" and "on2" 
distributions shown in Figures ^| and |^ have been nor- 
malized to match the area of the "off" region distri- 
butions. Areas around bright saturated stars (where 
the detection and photometry of faint stars is unre- 
liable) are excluded from the analysis. Table lists 
the usable area of the y-band image, properly ac- 



counting for the masked (excluded) regions around 
saturated stars. These masked regions are slightly 
different in the different bands for a given cirrus field; 
since saturation tends to be most severe at the longest 
wavelengths, the R- and /-band images have typi- 
cally smaller usable areas than the [/-band images. 
The saturated stars themselves are included in the 
star count analysis by adding the number of satu- 
rated stars in each field and each band to the bright 
end of the cumulative distributions. The total star 
count ranges from ~ 24,000 in U to ~ 51,000 in V 
for RCrA, and from - 2, 000 in U to ~ 7, 000 in V 
for Paley 1, the fields with the highest and lowest star 
densities, respectively, in our sample. 

A typical cumulative star count relation roughly 
resembles a power law, but has slight curvature such 
that the slope, S\ = dlog[iV(< m\))/dm\, gets shal- 
lower towards fainter apparent magnitudes (Figs. |^ 
and |^) . The relations are uncertain at the bright end 
due to the small number of stars; the abrupt flatten- 
ing at the faint end marks the onset of incomplete- 
ness. There is a range of slopes from field to field 
and a characteristic trend in S\ versus U BVRI, with 
the slopes steepening towards both shorter and longer 
wavelengths relative to the V band. This is likely the 
result of two competing effects. On the one hand, in- 
trinsically faint (lower main sequence) stars tend to 
have redder colors than intrinsically luminous stars 
(turnoff/upper main sequence), causing the 7-band 
slope to be steeper than the V-band slope. On the 
other hand, the "depth" of the sample, in terms of 
the typical distance of stars at the magnitude limit in 
each band, decreases from V towards U. The [/-band 
slope is close to the Euclidean value of 0.6 expected 
for a constant star density, while the mean distance 
of the F-band sample is large enough for the fall-off 
in Galactic star density along the line of sight to be 
noticeable. The star count slope is relevant to the 
discussion of the effects of dumpiness (Sec. 3.2). 

It is evident from Figures || and || that dust extinc- 
tion in the intervening cirrus cloud causes each on- 
cloud star count relation to be shifted towards fainter 
magnitudes (i.e., to the right) with respect to the cor- 
responding off-cloud relation. The magnitude of these 
shifts decreases from U to / as expected on the basis 
of the normal interstellar dust extinction law. More- 
over, the shift between the "on2" (the region with 
highest 100 (im flux) and "off" star count curves is 
consistently larger than the shift between "onl" and 
"off" curves across the UBVRI bands, at least in the 
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Table 2 
Properties of Cirrus Clouds 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Name 


(-FlOOjurrJbg 




Area off 


Area on i 


Area on2 


T 

on 




(MJysr- 1 ) 


(MJysr" 1 ) 


(deg 2 ) 


(deg 2 ) 


(deg 2 ) 


(K) 


RCrA 


9.3 


0.6 


0.61 


0.21 


0.15 


17.8 


PV1 


5.6 


0.3 


0.85 


0.43 


0.49 


17.5 


Paley 1 


0.19 


0.05 


0.47 


0.09 


0.32 


17.7 


Paley3 


1.0 


0.1 


0.83 


0.18 


0.16 


18.0 



Note.— 

Col. 2: Mean background 100 /im flux in regions without visible cirrus. 

Col. 3: Pixel-to-pixel r.m.s. variation in background 100 /tm flux. 

Col. 4-6: Area of "off", "onl", and "on2" regions (V-band), which refer 
to regions with fluxes that are < 5<r, 5er-llcr, and > IIct above the 100 
background. Note, these areas are slightly different across U BVRI because of 
differences in the area lost to saturated stars. 

Col. 7: Mean dust temperature in on-cloud region from DIRBE far infrared 
measurements (Schlegel et al. 1998). 
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two fields (RCrA and PV 1) for which the shifts are 
measured reliably for both "on2" and "onl" regions. 

4.2. Star Counts-Based Extinction Curves 

The magnitude shift between Wolf diagrams on 
and off an obscuring cloud is a direct measure of the 
effective dust extinction, A = ounts ; a t least at faint ap- 
parent magnitudes where the contribution of low lu- 
minosity foreground stars is unimportant and the on 
and off curv es run parallel to each other. As explained 
in Sec. 6.2 below, the y^ 01 " 1 ' 3 value derived from the 



magnitude shift in a Wolf diagram yields a measure 
of the effective "extinction" that depends on the de- 
gree of dumpiness of the dust distribution (even for 
a given mean dust optical depth and a given set of 
grain optical properties). 

Wolf diagram shifts are plotted in Figure | as a 
function of apparent magnitude in the UBVRI bands 
for all four cirrus cloud fields. The median value of 
^counts j g computed over an intermediate range of ap- 
parent magnitudes, avoiding excessive Poisson noise 
at the bright end and the divergence of the Wolf di- 
agram shift at the faint end caused by the flattening 
of the cumulative star counts due to incompleteness. 
In the cases of the "onl" and "on2" regions of the 
RCrA and PV 1 fields, the Wolf diagram shifts are 
sufficiently large for the extinction to be measurable 
in all five bands. For Paley 1 and Paley 3 on the other 
hand, the extinction is lower and the measurement 
errors larger (see discussion of errors below) so that a 
reliable measurement is only possible at the shortest 
wavelengths and that too for only the "on2" regions 
(U band for Paley 1, U and B bands for Paley 3); for 
these two fields the Ay and A B values are scaled to 
Ay assuming a specific form of the extinction law, ei- 
ther a standard law with Ry = 3.1 or an Ry — 1-7 
law (CCM). 

Figure |^ provides an illustration of the typical er- 
rors in the determination of the Wolf diagram appar- 
ent magnitude shift ^4™ unts . The errors in ^4c° unts are 
derived from the statistical uncertainties in the star 
counts, projected onto the magnitude axis (abscissa) 
using the slope of the log[A^(< m\))-vs-m\ relation 
(roughly a power law over the apparent magnitude 
range of interest), combined with the photometric 
uncertainties. The errors are of the order of 0.2 - 
0.3 mag for RCrA, the field with the largest number 
of stars. Poisson error dominates throughout most 
of the plotted apparent magnitude range; the contri- 
bution of the photometric error becomes noticeable 



beyond V ~ 17-18 in RCrA (somewhat fainter in 
Paley 3). Increasing photometric error, in conjunc- 
tion with the flattening of the cumulative star count 
relations, causes the flaring of the founts error bands 
at the faint end. 

The effective extinction curves for the "onl" and 
"on2" regions of RCrA and PV 1 derived from star 
counts are shown as open squares and triangles, re- 
spectively, in Figure ^. The yl™ unts estimates are un- 
certain for the Paley 1 and Paley 3 fields (all but the 
"on2" estimates at the shortest wavelengths) so ex- 
tinction curves are not presented for these two clouds. 
The lines are analytic approximations to the extinc- 
tion law provided by CCM for different Ry values. 
The RCrA and PV 1 extinction curves are reasonable 
approximations to the standard extinction law with 
Ry =3.1 (solid line). 

Looking more closely, the PV 1 "on2" extinction 
curve rises more steeply towards short wavelengths (U 
band) than the standard extinction law, suggesting a 
lower value of Ry. A minimum of the x 2 statistic is 
obtained for Ry = 1.7 (dashed line). For RCrA on 
the other hand, the increase of founts with decreas- 
ing wavelength is somewhat more gradual than for the 
standard extinction curve, with ^4™ unts being partic- 
ularly discrepant. While this would appear to favor 
high Ry values in that the \ 2 statistic decreases with 
increasing Ry, there is no clear \ 2 minimum over the 
plausible range of parameter space (Ry < 7). Alter- 
natively, a good fit to the RCrA founts measuremen ts 
is obtained by assuming a non-uniform dust density 
(dotted line; Sec. |6 . 2[ ) with Ry = 3.1 or even with 
Ry = 1.7, the lower Ry value being preferable in light 
of the color excess data for this cloud (see Sec. 5.1). 



We will return to a discussion of the counts-based ex- 
tinction curves in Sec. 



3.2 



The weighted average of the founts measurements 
in UBVRI is computed for RCrA and PV1, after 
the measurement in each band is scaled to the equiv- 
alent Ay estimate using the analytic form of the ex- 
tinction law for a specific Ry (Table || and Fig. |l3| ). 
While the choice of Ry affects the quality of the fit to 
the five ^4™ unts measurements, the weighted average 
^counts j g p rac tically independent of Ry , and tends to 
be similar to the measured value of Amounts in the V 
band. In other words, changing Ry results in a steep 
or shallow extinction curve, but the best-fit curves in- 
tersect near the middle of the range of available data 
points which happens to be the V band. 



9 



4.3. Star Count Models 

Wolf diagrams have traditionally been used to de- 
termine both the extinction through and the distance 
to dust clouds. Distance determination, in particular, 
requires knowledge of the stellar luminosity function 
and the spatial density distribution of stars along the 
line of sight. Based on this information, theoretical 
Wolf diagrams can be constructed for various values 
of the dust extinction and distance of an obscuring 
cloud. These diagrams can then be directly compared 
to the star count data. 

Theoretical Wolf diagrams are generated for the 
four cirrus clouds in this study using the IASG Galaxy 
Model program (Bahcall, Casertano, & Ratnatunga 
1987; Ratnatunga, Bahcall, & Casertano 1989; Caser- 
tano, Ratnatunga, & Bahcall 1989), which is based on 
the Bahcall-Soneira star count model of the Galaxy 
(Bahcall & Soneira 1980; Bahcall 1986). The IASG 
program yields the total projected density of stars 
and their distribution as a function of apparent V 
magnitude, B — V color, and distance modulus, for 
any specified direction. A Monte-Carlo realization of 
these distributions is carried out in order to construct 
a stellar sample. Next, U — B colors are computed us- 
ing the Bertelli et al. (1994) stellar evolution models, 
averaging over all available ages per metallicity, for 
four different metallicities, Z = 0.0004, Z = 0.008, 
Z = 0.05 and Z = 0.02. The product is a list of 
stars with UBV magnitudes and distance moduli for 
each metallicity from which cumulative star counts 
are computed appropriate for the off-cloud region. 
The dust cloud is assumed to be (geometrically) thin; 
it produces an apparent offset of A\ for all stars lo- 
cated at distances greater than that assigned to the 
dust cloud. The resulting list is used to construct 
on-cloud model star count relations. As the model 
results are similar for the four different metallicities, 
only curves computed with Z = 0.05 are shown. 

The slopes of the model cumulative star count rela- 
tions, plotted as dotted blue lines in Figures ^ and ^, 
agree well with the observed slopes, with the possi- 
ble exception of P V 1 . An offset has been applied to 
the apparent magnitude scale in the Bahcall-Soneira 
model in order to match the observed counts in the 
"off" region: +0.2, +1.5, -0.2, and -0.5 mag for 
RCrA, PV1, Paley3 and Paleyl, respectively. Pos- 
itive corrections to the apparent magnitude may in- 
dicate that the overall extinction integrated over the 
line of sight is not quite zero even in the off-cloud re- 



gions. There may, however, also be small systematic 
errors in the model (e.g., in the apparent magnitude 
scale, or in the normalization of the stellar density or 
luminosity function) that cause it to deviate from the 
observed counts. Note that since the curves (model 
and observed counts) are good approximations to sin- 
gle power laws, a horizontal shift is roughly equivalent 
to a vertical shift [in the log(iV) direction] . 

4.4. Cloud Distance Determination 

Model on-cloud and off-cloud Wolf diagrams are 
constructed for the RCrA and PV 1 fields using IASG 
star count predictions coupled with a geometrically 
thin dust slab located at various line-of-sight dis- 
tances. These model curves are shown as dotted blue 
lines in the upper panel of Figures ^ and || for a cloud 
distance of 550 pc and an optical depth of Ajj = 0.9 
(matching the "on2" measurements of A^ unts in the 
two clouds), which scales to Ay — 0.6 on the basis of 
a standard extinction law. The on-cloud model curves 
for RCrA shown in the center and bottom panels of 
Figure are for cloud distances of 400 pc and 800 pc, 
respectively, for optical depths of Ajj = 0.9 ("on2"; 
left) and A v = 0.3 ("onl"; right). 

Only RCrA and PV1, the fields with the low- 
est latitiude and largest area, respectively, contain 
enough bright stars (U ~ 13 - 15; see discussion be- 
low) and high enough extinction for us to obtain re- 
liable constraints on the distance using this method. 
Figure ^ compares the RCrA "on2" and "onl" Wolf 
diagrams (upper left and upper right panels, respec- 
tively) with their corresponding IASG model-based 
Wolf diagrams for cloud distances of d = 0.4 kpc (cen- 
ter panel) and d — 0.8 kpc (bottom panel), bracketing 
the likely distance of the RCrA cloud. The best es- 
timate of the RCrA cloud distance is in the range 
550- 750 pc. For PV1, it is only possible to place 
an upper limit to the distance of about 1 kpc. The 
distance determination is uncertain for PV 1 due to 
its low stellar surface density: the Poisson error in 
this field is 40% larger than in RCrA despite its 3.3 x 
larger area. 

The IASG model Wolf diagrams shown in the cen- 
ter and bottom panels of Figure |?] are based on a very 
large number of synthetic stars in order to minimize 
the effect of Poisson error. A realistic measure of the 
dispersion in cumulative star count relations is ob- 
tained by making several Monte Carlo realizations of 
the on-cloud model, with each realization containing 
the number of stars appropriate for the actual area 
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of the RCrA on-cloud ("on2" or "onl") region being 
simulated. The dispersion among the realizations of 
a given model can be used to judge the significance 
with which the model fits or is ruled out by the data; 
the ±lcr error bars in the center left panel of Fig. ^ in- 
dicate the dispersion for the RCrA "on2" d = 0.4 kpc 
model. Averaging over the U = 12-15 mag range, 
66% of the d — 0.4 kpc model realizations lie be- 
low the RCrA "on2" star count curve, while 82% of 
the d — 0.8 kpc model realizations lie above it (both 
models adopt Ajj = 0.9). Comparing the Ajj = 0.3 
model to the RCrA "onl" star count curve, 54% of 
the d = 0.4 kpc realizations lie below the data while 
90% of the d — 0.8 kpc realizations lie above the data. 

The differences between the on-cloud and off-cloud 
model cumulative star count curves can be under- 
stood in simple terms. The transition between majority- 
foreground and majority-background in the on-cloud 
star count relation occurs at an apparent magnitude 
equal to that of a "typical" star at the distance of 
the cloud. While stars have a wide range of absolute 
magnitudes, the characteristic absolute magnitude is 
defined by an abrupt rise near the bright end of the 
stellar luminosity function. The rise is at the main 
sequence turnoff, M^? STO w +3.5, for any old stellar 
population, and is particularly steep in the U band 
(M^ 1 STO w +4.0). Thus a cirrus cloud at a distance 
of 0.4 kpc (0.8 kpc), corresponding to a distance mod- 
ulus of (m — M)o = 8.0 mag (9.5 mag), is expected to 
have a transition region at U ~ 12.0 (U ~ 13.5), in 
agreement with the model curves shown in the center 
(bottom) panel of Figure [?]. 

5. Effect of Reddening on Stellar Colors 

In this section, the distribution of stellar colors is 
compared between on- and off-cloud regions for the 
purpose of exploring the reddening (selective extinc- 
tion) caused by dust. In contrast to the star count 
analysis described in the previous section, the color 
analysis described in this section makes no use of the 
number or number density of stars in the various re- 
gions, and only uses information about the normal- 
ized stellar color distribution in each region ("off", 
"onl", and "on2"). As discussed in Sec. ^| below, star 
counts and color distributions provide complementary 
measures of the optical depth of a cirrus cloud. 



5.1. Cumulative Color Distributions and Color 
Excess Measurements 

A direct way to determine the amount of reddening 
caused by the cirrus cloud is to measure the difference 
in the stellar colors between on- and off-cloud regions. 
This is done by constructing normalized cumulative 
color histograms in U — V, B — V, R — V, and I — V 
for the "off" , "onl" , and "on2" regions of the cirrus 
cloud fields. The median horizontal (i.e., color) dif- 
ference between on-cloud and off-cloud histograms is 
defined to be the color excess, E(m\ — V) = A\ — Ay, 
for A = U BRI. Sample cumulative U — V and I — V 
color histograms are displayed in Figure || (solid and 
dashed lines, respectively) for the "on2" and "off" re- 
gions (thin and bold lines, respectively) of RCrA (up- 
per panel) and Paley3 (lower panel), with the hor- 
izontal line segment indicating the median value of 
the offset in each case. The color excess is well deter- 
mined for most fields and colors, with the exception 
of long wavelength measurements for the low optical 
depth clouds, Paleyl and Paley3 [as illustrated by 
the Paley3 E(I - V) data in Fig. §]. 

If a small fraction of the area of the "onl" or "on2" 
region of the cloud happens to have a much larger op- 
tical depth than the rest, it is effectively ignored by 
the use of the median in the color offset computa- 
tion; the significance of this aspect of our measure- 
ment technique is discussed further in Sec. 6.2. The 



median offset is also insensitive to dilution by unred- 
dened, bright, foreground stars, especially as these are 
greatly outnumbered by background stars. 

The median color excess, E(m\ — V) for m\ = 
UBRI, for the "onl" and "on2" regions of the four cir- 
rus clouds is plotted as open squares and triangles in 
Figure ^|. The dashed and solid lines show the CCM 
analytic approximation for the Galactic interstellar 
reddening law for Ry = 3.1 for "onl" and "on2", re- 
spectively, while the dotted line shows an Ry = 1.7 
reddening law for "on2" . Each analytic curve is scaled 
by the appropriate weighted mean value of the visual 
optical depth, A™ ccss , derived from the four color ex- 
cess measurements (listed in Table |^ and shown in 
Fig. |l^) . As expected, ^4^ cess is sensitive to the choice 
of Ry; this is in contrast to the ^4™ unts estimate de- 
rived from the founts measurements (Sec. |4.2| ). 

The measured color excesses tend to increase more 
steeply with increasing 1/A (V — > U) than the best- 
fit standard Ry = 3.1 reddening law, and vary less 
strongly with A from V — ► / than the standard red- 
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dening law. In all cases, the Ry — 1.7 reddening law 
shown in Figure || appears to be a better match to the 
shape of the observed E(m\ — V) versus 1/A than the 
standard reddening law. Applying a x 2 test to the 
RCrA and PV 1 color excess data, a minimum of the 
X 2 statistic, albeit a relatively shallow one, is found 
for Ry values that are substantially lower than the 
standard Galactic value of 3.1. The Paleyl and Pa- 
ley 3 x 2 behavior also favors a low value of Ry (<; 2) 
but there is no clear minimum as the Poisson error in 
these fields is larger than for RCrA or P V 1 . 

Note, the x 2 test is only used in a relative way here. 
The x 2 values are not used to decide if the CCM an- 
alytic form of the reddening law with a specific value 
of Ry is a good fit to the data in an absolute sense; 
the exact value of Xmin ^ s n °t particularly meaningful 
because: (1) the E(m\ — V) error estimates used in 
the x 2 test are only approximate representations of 
the true (likely non-Gaussian) error distribution; and 
(2) the Ry -parametrized analytic formula used in the 
analysis is known to deviate from the measured red- 
dening law in certain lines of sight (Cardelli & Clayton 
1991). 

5.2. Color-Color Diagrams 

A second way to determine the cloud optical depth 
is through the use of color-color diagrams. This 
method is similar in principle to the color excess 
method described above except that a pair of stel- 
lar colors [e.g., the (U — B, B — V) pair or the 
(B — V, V — I) pair] is analyzed jointly. As the in- 
trinsic spread in stellar V — R and R — I colors is 
not very broad relative to the photometric errors in 
these quantities, the i?-band data are excluded from 
this analysis and the longer V — I color baseline is 
used instead for a more accurate measurement of the 
reddening. 

Figure |l^ shows the U — B versus B — V and 
V — I versus B — V color-color diagrams (upper and 
lower panels, respectively) for RCrA and PV1 (left 
and right panels, respectively). The Paleyl and Pa- 
ley 3 fields contain too few stars and the reddening 
in these clouds is too low for an analysis of their 
color-color distributions. The black, red, and green 
points indicate stars in the "off" , "onl" , and "on2" 
regions, respectively, plotted in this order. The blue 
lines show theoretical loci, computed from the Bertelli 
et al. (1994) stellar evolution models by averaging 
all available ages per metallicity, for three metallic- 
ities: Z = 0.0004, Z = 0.008, and Z = 0.05 (top 



to bottom). The large spread of points is due to a 
combination of photometric errors, intrinsic range of 
stellar colors (spread in metallicity and age), sample 
contamination by background galaxies, spurious de- 
tections/matches, etc.. 

The arrows in Figure ^| indicate reddening vec- 
tors (CCM) corresponding to Ay = 2 mag for a stan- 
dard extinction law with Ry =3.1 (red) and for an 
Rv — 1.7 extinction law (black). The latter extinc- 
tion law is a reasonable approximation to the available 
star count and color excess data for the four clouds. 
The progressive shift along the reddening vector of 
the "onl" and "on2" stellar color-color distributions 
with respect to the "off" distribution is clearly visi- 
ble in RCrA, and is somewhat less obvious in PV1. 
The principal direction along which the stellar distri- 
bution is extended is almost parallel to the reddening 
vector, especially in the V — I vs B — V plot; even 
so, an accurate statistical measurement of the redden- 
ing is rendered possible by the use of a well-matched 
off-cloud region and the large numbers of stars in the 
sample (which enables accurate determination of the 
centroid of the stellar distribution). 

The amount of reddening in the color-color dia- 
grams is estimated by measuring the shifts between 
on- and off-cloud stellar distributions projected onto 
the direction of the reddening vector. The sample of 
stars is divided into a series of contiguous slices paral- 
lel to the reddening vector each of width 0.1 mag. The 
median offset along the length of the slice is computed 
between the "onl"/"on2" and "off" distributions for 
each slice. The final visual optical depth estimate is 
a weighted average of the median shifts derived from 
the different slices. These optical depth estimates are 
denoted Ay BV and A$ VI for the (U — B, B — V) 
and (V — I, B — V) plots, respectively. The offsets 
(in Ay units) for the individual slices of the RCrA 
U — B vs B — V diagram are illustrated in Figure 11 
for "on2" (solid line) and "onl" (dashed line). The 
error bars indicate the combined standard error in the 
mean centroids of the on- and off-cloud stellar distri- 
butions. The weighted mean of the individual slice 
offsets is indicated by the horizontal dot-dashed line. 

A test is carried out in order to assess the ef- 
fect of incompleteness at the faint end of the sample, 
or, more specifically, the effect of differential incom- 
pleteness between on- and off-cloud samples which 
may introduce a bias in the color offset measurement. 
The "off" sample is artifically reddened by an arbi- 
trary amount corresponding to an optical depth of 
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Ay m . This simulated reddened sample is then folded 
through the completeness function in the three bands 
of interest (i.e., UBV or BVI). The offset of the on- 
cloud sample is measured with respect to the simu- 
lated reddened and truncated sample. The process is 
iterated while varying Ay m until there is a null offset 
between simulated and on-cloud samples. The value 
of A v m that produces a null offset is found to be vir- 
tually identical to the directly measured value of the 
visual optical depth (A V BV or A V VI ), so it may be 
concluded that differential incompleteness effects are 
unimportant in the determination of Ay from color- 
color diagrams. 

The calculation of A V BV and A B VI for the "onl" 
and "on2" regions of RCrA and PV 1 is done two ways, 
assuming a standard Ry = 3.1 reddening law or 
the steeper Ry = 1.7 reddening law. The result- 
ing Ay estimates are listed in Table || and shown in 
Figure [IJ|. The second set of Ay values are under- 
standably lower than the first: a lower Ry produces 
a steeper extinction law in the blue portion of the 
spectrum (i.e., stronger variation with A), so that a 
fixed (observed) color offset corresponds to a smaller 
Ay. A secondary effect is that the directions of the 
Ry = 3.1 and Ry = 1.7 reddening vectors are slightly 
different in the (V — I, B — V) plane and this leads to 
different projections of the stellar color-color offsets. 

In contrast to the star count and color excess anal- 



ysis described above (Sec. |4.2| and 5.1), a specific form 
of the reddening law had to be assumed a priori in de- 
riving Ay from color-color diagrams. Given the sim- 
ilar orientations of the Ry = 3.1 and Ry = 1.7 red- 
dening vectors (directions almost degenerate in U — B 
vs B — V plot), it is impossible to determine both Ry 
and Ay independently from a single color-color di- 
agram, even for a field containing as many stars as 
RCrA. It would require a very large, uncontaminated 
sample of stars and accurate photometry to exploit 
the subtle dependence of reddening vector direction 
on Ry in the (V — I, B — V) plane. Alternatively, the 
combination of all four available colors ca n di scrimi- 
nate between different Ry values (cf. Sec. |5.l|) . 

6. Discussion 

6.1. Reddening Law in Diffuse Cirrus Clouds 

Studies of the interstellar extinction law along vari- 
ous lines of sight in the Galaxy show some evidence for 
a dependence of Ry on environment. Dense regions 
such as optically thick molecular clouds, of which 



the p Oph dark cloud is a typical example (Martin 
& Whittet 1990), tend to have relatively shallow ex- 
tinction curves indicating Ry values larger than the 
global average of 3.1, while regions of low density dis- 
play extinction curves that rise steeply towards short 
wavelengths indicating lower-than-average Ry values 
(cf. Mathis 1990), although some counter examples 
do exist. The fact that all four diffuse clouds in our 
sample have reddening laws that favor lower Ry val- 
ues is in keeping with this trend. Also consistent with 
the findings of this paper is the observed steepening 
of the far-ultraviolet extinction curve with increasing 
distance from the Galactic plane (Kiszkurno-Koziej 
& Lequeux 1987). On the face of it, the RCrA data 
appear to yield conflicting results: a higher-than- 
average Ry value fits the star counts-based extinc- 
tion curve (Sec. f|) while a lower-than-average value 
fits the color exc ess measurements. As explained in 
Sec. 4.2 and |6.2|, however, all of the RCrA data are 



adequately explained by Ry <J 2 if one invokes a non- 
uniform dust density. 

The Ry values for the four cirrus clouds in this 
study are unusually low, although comparable values 
have been observed in a few instances. A study of 
about 200 OB stars by Whittet & van Breda (1980) 
found a minimum value of Ry = 2.4, while the sub- 
sequent survey by Fitzpatrick & Massa (1990) found 
a minimum value of Ry = 2.6 for the line of sight 
towards HD 204827. More recently, Larson, Whittet, 
& Gough (1996) reported a value of Ry = 2.1 for the 
high-latitude cloud DBB 80 based on their study of 
the early type star HD 210121. 

The observed variations in the shape of the inter- 
stellar extinction law from one line of sight to another 
are predominantly at blue and ultraviolet wavelengths 
while the red and infrared portion remains invariant 
(Martin & Whittet 1990). Variations in Ry, the ra- 
tio of selective blue-minus-visual extinction to total 
extinction, are generally thought to be a result of 
variations in the size distribution of dust grains: the 
total visual optical depth is determined by the abun- 
dance of large grains, while the extinction in the blue 
and ultraviolet portion is dominated by smaller grains 
(cf. CCM; Larson et al. 1996). Thus, a steep extinc- 
tion curve is caused by a higher than average propor- 
tion of small size particles in a dust cloud. An abun- 
dance of small particles may be explained in terms of 
the cirrus being in a UV-poor and weak radiation en- 
vironment in which grain destruction by photon heat- 
ing is inefficient. An alternative explanation is that 
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the growth, by coagulation, of larger size particles is 
inhibited by the low overall density of the cloud (Lar- 
son et al. 1996). 

6.2. Effect of a Clumpy Medium 

In this section, we investigate the effect of a clumpy 
cirrus cloud on the counts and colors of background 
stars. The technique of using normalized color his- 
tograms to measure the reddening caused by cirrus 
(Sec. ||) yields the average optical depth for the lines 
of sight towards the specific set of stars used in the 
color analysis. If there is substantial non-uniformity 
in the dust column density within a given cloud re- 
gion ("onl" or "on2"), so much so that a highly- 
reddened subset of background stars drops out of 
the matched U BVRI sample, the color-based method 
will underestimate the true area-weighted average op- 
tical depth. Wolf diagrams, on the other hand, are 
sensitive (in a statistical sense) to stars that have 
dropped out of the on-cloud sample (Sec. |J). In addi- 
tion to the "horizontal" shift to fainter apparent mag- 
nitudes of the on-cloud cumulative star count curve 
due to extinction of stars in the sample, any highly- 
attenuated undetected subset of stars will cause a 
"vertical" downward shift in the curve, which is equiv- 
alent to a shift to the right since the star count rela- 
tion closely resembles a power law. 

The characterization of a clumpy medium is, in 
general, complicated and beyond the scope of this pa- 
per. Instead, the essential features of a non-uniform 
dust cloud are illustrated with the help of a simple 
two-phase model: a fractional area / of the absorbing 
material is assumed to have a uniform and relatively 
low optical depth of cct while the remaining frac- 
tion (1 — /) is covered by dense clumps with optical 
depth ^4y Ump . Exact calculations are carried out for 
the particular case: Ay Ump — ► oo. Since the Wolf 
diagram shift is statistical in nature, the following 
calculations involve no assumptions about the mor- 
phology of the low and high extinction regions. How- 
ever, considering that a large contiguous area with 
Ay — > oo would be conspicuous in the CCD images, 
a low extinction sheet with numerous small dense 
clumps seems a more likely scenario. 

With the above assumptions, and given the specifics 
of the measurement method used (Sec. ||), the color 
excess histograms and color-color diagrams measure 
the extinction in the smooth, low optical depth part 
of the cloud, A x hcct . The Wolf diagram shift is greater 



than ^4 A hcot but its exact value depends on the shape 
of the cumulative star count function. For the pur- 
poses of this calculation, the off-cloud number count 
relation is assumed to be a power law with slope 
S\: N oS {< m x ) = C x 10 s ^ m \ where C x is a nor- 
malization constant. The cumulative star count func- 
tion in the on-cloud ("onl" or "on2") region is then 
a power law with the same slope: N on (< m\) = 



C 



Thus. 



f ioSaK-aD + (l - /)io*a< 

the Wolf diagram shift, (A™ unts ) modc i, is directly pro- 
portional to the "vertical" shift between the on- and 
off-cloud relations at a fixed apparent magnitude m\. 



1 



(A counts )modci = - [logoff) - log(AU)] (3) 



log 



/10 



x +(1-/)10 



-SxA": 



(4) 



Under the assumption that A^ ump » 1 over the 
wavelength range of interest, this simplifies to: 



Ddcl 



log(/) 
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sheet 



(5) 



where the first term accounts for the undetected frac- 
tion of background stars while the second term is the 
extinction caused by the smooth sheet on the rest of 
the stars. The I/-band version of the above equation 
may be rewritten as: 



log(/) 



S v [Af eet 



{ a counts \ 
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(6) 



Normalizing the counts-based "effective" extinction 
(Eq. ||) by the y-band value and substituting for / 
using Eq. [|, it follows that: 



^counts 



^counts 



jdcl 
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A sheet 
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( A counts^ 

^ model 



(7) 



where r\ = (A\/Ay) depends only on the shape of 
the extinction law. The UBVRI star count slopes 



vary somewhat from field to field (see Sec. 4.1 for a 
discussion of the general trends in 5a). 

A one-parameter family of normalized counts-based 
"effective" extinction curves, computed by varying 
the ratio A^ cct /(A^ ounts ) mod ci for a standard Ry = 
3.1 extinction law (CCM) and using the observed 
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star count slopes S\ in the RCrA field, is shown in 
Figure 12. An increasing degree of non-uniformity 



(decreasing /) results in a flattening of the counts- 
based extinction curve. The effect of dumpiness is 
qualitatively similar for other extinction laws as well 
(e.g., Ry = 1-7). Although Eq. [?] depends only on the 
ratio A^ eet /A^ unts and not on the value of A c J nnts it- 
self, the derived value of the area covering fraction / 
does depend on Ay unts (Eq. Moreover, the nor- 
malized counts-based extinction curves display a sim- 
ilar behavior with increasing dumpiness even if the 
^ciump _^ ^ assumption is relaxed, although their 
exact shapes in that case depend somewhat on the 
value of ^^ ounts . 

There is an apparent discrepancy in the case of 
RCrA. The Wolf diagram shifts for this cloud, A c ° unts , 
increase less steeply with decreasing wavelength than 
the standard Ry = 3. 1 e xtinction law, indicating a 
larger Rv value (Sec. 4.2). On the other hand, the 
color excesses for the same cloud are a steep function 
of wavelength and are adequately described by the 
Rv = 1.7 CCM reddening law (Sec. |t|). This dis- 
crepancy can be resolved if the dust column density in 
RCrA is non-uniform. For example, the founts mea- 
surements can be fit by a combination of Ry — 1.7 
and a smooth sheet fraction of / = 0.8, which corre- 
sponds to A^ cct /A™ unts = 0.5 and Af cct = 0.3 (this 
translates to / = 0.7 for A™ unts = 1 in Fig. 0). The 
combinations (Rv = 3.1, / = 0.9) and (Rv ~ 5, / = 
1) produce ^4= ounts vs A functions that are practically 
identical in shape to the (Ry = 1.7, / = 0.8) case 
and they all fit the RCrA measurements equally well 
(see dotted line in Fig. |(]). Considering the low Ry 
value implied by the RCrA color excess data, the first 
solution seems the most plausible, in spite of the de- 
generacy of the fits to the founts ^ata. 

6.3. Comparing Different Measures of the 
Dust Optical Depth 

The preceding sections outline a variety of methods 
used to measure the dust extinction optical depth of 
a few representative diffuse interstellar cirrus clouds: 
empirical scaling of the 100 fxm. flux after correction 
for temperature variations (Sec. ||), analysis of the 
background star count density as a function of ap- 
parent magnitude (Se c. ||), reddening of stellar color 
distributions (Sec. 5.1), and shifts of the stellar color- 
color distribution (Sec. |5.2| ). The different methods 
measure slightly different quantities (as quantified in 
Sec. 5.2). Nevertheless, the methods yield general 



consistency for a given cloud. The full set of Ay es- 
timates, based on both the standard Ry = 3.1 value 
and the preferred Ry = 1.7 value, are listed in Table || 
and plotted in Figure |l3[ 

The agreement between different Ay estimates for 
the PV 1 cloud is especially good if one adopts Ry <; 2 
in deriving Ay from the color excess and color-color 
data (filled symbols in the second panel of Fig. [l3]) — 
the observed trends of A™""* 8 and E(m\ — V) versus A 
for this cloud are, in any case, indicative of a low Ry. 
The RCrA "on2" region Ay estimates, particularly 
those based on Ry — 1.7 (filled triangles in panel 1 
of Fig 



13), display the largest discrepancy: 



A counts 
A V 

is twice as large as estimates based on the color ex- 
cess and color-color distribution and is significantly 
smaller than the 100 urn flux-based estimate, Ay 1 0Aim . 
As discussed in Sec. 3^ above, this is best explained 
in the context of the simple two-phase model of a 
clumpy cirrus medium. The dotted line in panel 1 
of Figure [l3| indicates a model designed to match 
the RCrA "on2" observations: Ry = 1.7, / = 0.8, 

^sheet = o 3j and ^counts _ q g egtimates 

based on color excess and color-color data, Ay* 0085 , 
A V BV , and A V VI , correspond to Af cct . The counts- 
based estimate, Ay unts , is boosted by the presence of 
high density clumps. The far infrared flux-based esti- 
mate, Ay^ is expected to be larger than both the 
color- and counts-based estimates and is nominally 
a weighted average of A^ oct and A v ump . However, 
since the characteristic size of the dense clumps is 
likely to be smaller than the coarse (~ 1°) angular 
resolution of the far infrared data and since the dust 
in the dense clumps is likely to be colder than the dust 
in the low density sheet (because of self-shielding), the 
flux-weighted average temperature derived by SFD is 
likely to be biased in favor of the temperature of the 
warmer dust in the sheet. The Ay estimate derived 
from the scaling of i^iooLm (Sec. |J) is then lower than 
the "true" area-weighted average optical depth. The 
lower and upper branches of the dotted line in panel 1 
of Figure Il3| show the area-weighted average Ay for 



.4' 



clump 
V 



2.4 and 3.5, respectively; these bracket the 
on2". The actual opti- 
cal depth of the dense clumps could be significantly 
higher, so it is plausible that stars behind such clumps 
drop out of the matched UBVRI sample — this is 
equivalent to settin g A v nmp — > oo in the Wolf dia- 
gram analysis (Sec. |6.2| ). 

The ratio of visual optical depth to far infrared 
brightness, Ay / -Fioofuni is in the range 0.06-0.08 for 



measured ^4^°^ m f or RCrA 
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Table 3 
Extinction Estimates 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Name 


Rv 


Region 




A counts 


4 excess 


A V BV 


Ay VI 



RCrA 


3.1 


onl 


0.25 ± 0.02 


0.26 ± 0.04 a 


0.10 ±0.02 


0.18 ±0.01 


0.05 ±0.01 




3.1 


on2 


0.77 ±0.04 


0.61 ± 0.04 a 


0.46 ± 0.04 


0.45 ± 0.01 


0.36 ±0.01 




1.7 


onl 




0.26 ± 0.03 a 


0.07 ±0.01 


0.11 ±0.01 


0.06 ±0.01 




1.7 


on2 




0.61 ± 0.04 a 


0.27 ±0.02 


0.26 ± 0.01 


0.25 ±0.01 


PV1 


3.1 


onl 


0.16 ±0.01 


0.24 ±0.06 


0.33 ±0.18 


0.36 ± 0.04 


0.28 ±0.04 




3.1 


on2 


0.35 ± 0.01 


0.46 ± 0.07 


0.65 ±0.18 


0.55 ± 0.03 


0.52 ±0.04 




1.7 


onl 




0.21 ±0.06 


0.20 ±0.11 


0.16 ±0.02 


0.13 ±0.03 




1.7 


on2 




0.42 ± 0.06 


0.37 ±0.10 


0.31 ± 0.02 


0.32 ± 0.02 


Paley 1 


3.1 


onl 


0.05 ±0.01 




0.16 ±0.06 








3.1 


on2 


0.16 ±0.01 


0.12±0.23 b 


0.24 ±0.07 








1.7 


onl 






0.11 ±0.04 








1.7 


on2 




0.09±0.17 b 


0.16 ±0.04 






Paley 3 


3.1 


onl 


0.06 ±0.01 




0.08 ±0.08 








3.1 


on2 


0.14 ±0.01 


0.13±0.11 b 


0.17 ±0.08 








1.7 


onl 






0.06 ±0.05 








1.7 


on2 




0.09 ± 0.09 b 


0.10 ±0.04 







accounts j g determined by fitting the U BVRI Wolf diagram shifts to a clumpy cirrus cloud model: standard Rv =3.1 
extinction law with a fr actio n (1 — /) = 0.1 of the cloud area covered by optically thick dense clumps, or with Rv = 1.7 
and (1 - /) = 0.2 (Sec. 

b^counts ig scale d, using the standard R v = 3.1 extinction law or R v = 1.7, from A c ^ unts in the Paley 1 "on2" region 
or from y^° unts and A c £ unts in the Paley 3 "on2" region. The rest of the "on2" and all "onl" Al° unts values are poorly 
constrained. 



Note. — 

Col. 4: Derived from 100 /im flux (Sec. |). 

Derived from UBVRI star counts/Wolf di agra ms (Sec. 
Derived from UBVRI color excesses (Sec. 
Derived from UBV color-color diagram (Sec. 
Derived from BVI color-color diagram (Sec. 
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RCrA ("onl" and "on2"), Paleyl ("on2"), and Pa- 
ley3 ("on2"), and in the range 0.12-0.13 for PV1 
("onl" and "on2"). These results are in good agree- 
ment with other measurements of diffuse interstellar 
clouds: 0.05-0.1 in Stark's (1995) study of isolated, 
high-latitude clouds (including Paleyl and Paley3); 
0.05 in an all-sky survey by Boulanger & Perault 
(1988); and 0.04 in a cirrus cloud in the direction 
of the Polaris star (Zagury et al. 1999). 

The -Fioo^m values quoted above for the cirrus 
clouds in this study are the calibrated far infrared 
fluxes from SFD, uncorrected for dust temperature. 
Using the temperature-corrected flux, F^q" (the 
flux a cloud would have if it were at the global av- 
erage dust temperature of (T) = 18.0 K), lowers the 
A v /Ff^ m ratio to 0.04-0.06 for RCrA, Paleyl, and 
Paley3, and to 0.07-0.09 for PV1. These clouds 
are colder than the typical dust cloud in the Galaxy 
(see Table ^|), probably due to the absence of strong 
heating sources. Thus, the temperature correction 
causes F^fi™ to be greater than Fioo^m, and the 
Ay /i^xoSum to De correspondingly lower. The typ- 
ical large dust cloud complex in the Boulanger & 
Perault (1988) study is in a hotter radiation envi- 
ronment than the cirrus clouds in our sample and the 
Zagury et al. (1999) cloud is thought to be heated in 
part by the Polaris star; this may explain why their 
Av/FiQQpja measurements are at the low end of the 
range found in this study. If the SFD temperature 
estimates are taken literally, the observed difference 
between Ay / i*xoo/im for P V 1 versus the other clouds 
is not due to cloud-to-cloud variations in dust temper- 
ature, but instead reflects variations in dust optical 
properties and/or measurement error. 

6.4. Future Prospects: Sloan Digital Sky Sur- 
vey Data 

The techniques described in this paper are applied 
to four representative cirrus clouds. However, the ap- 
plication and the results derived from it are limited 
by several factors including: field of view of the ob- 
servations (~ 1 deg 2 ), photometric depth of the stel- 
lar sample (V5i m « V50 ~ 20), photometric accuracy 
(typically a ^ 0.1 mag), and the cirrus cloud sam- 
ple size (4). In the case of very thin cirrus, a large 
number of background stars and accurate photome- 
try are needed to make an accurate determination of 
the extinction. This can be accomplished by using a 
larger field of view and/or increasing the photometric 
depth. Determination of cloud distance can only be 



achieved with a substantial number of (bright) fore- 
ground stars, and this requires a wide field of view. 
A large sample of optically thin cirrus clouds must be 
studied before the conclusions about low Ry, small 
grains, and dumpiness can be extrapolated to the 
general diffuse ISM. 

The Sloan Digital Sky Survey (SDSS) data set 
should combine several of the desired properties out- 
lined above. The SDSS is an ongoing digital pho- 
tometric and spectroscopic survey covering 10 deg 2 
in the North Galactic cap, and a deeper imaging 
survey in the South Galactic hemisphere covering 
225 deg 2 . Accurate photometry is being carried out 
in five bands, v! , g 1 , r 1 , i' , and z' (spanning the entire 
near-ultraviolet to optical range), complete to limit- 
ing point source magnitudes of uiab = 22.3, 23.3, 
23.1, 22.3, and 20.8, respectively in the North Galac- 
tic cap. The southern portion of the survey should 
extend about 2 mag fainter than the northern por- 
tion. The entire photometric survey is expected to 
include over 10 8 stars and galaxies. 

Application of the photometric methods explored 
in this paper to the SDSS database should yield a re- 
liable and detailed extinction map, one that is based 
on direct measurement as opposed to scaling of far in- 
frared flux, along with excellent statistics for the mea- 
surement of cloud distances. This exercise will likely 
go hand in hand with the development of a new em- 
pirical star count model of the Galaxy, a refinement of 
the Bahcall-Soneira model. The SDSS northern sam- 
ple is expected to be about 2 mag deeper than the 
sample analyzed in this study, and this corresponds 
to about an order of magnitude increase in the surface 
density of background objects. The increased density 
of background tracers should translate to an increase 
in the angular "resolution" of the resulting extinction 
map, especially for the southern part of the SDSS sur- 
vey. The SDSS opens up the prospect of obtaining a 
three-dimensional map of cirrus clouds over a quar- 
ter of the sky and conducting an extensive study of 
diffuse interstellar dust clouds. 

7. Summary 

We develop techniques for investigating the opti- 
cal depth, constituent grain properties, structure, and 
distance of diffuse interstellar dust clouds. The study 
is based on UBVRI CCD photometry of several thou- 
sand stars in four > 1 deg 2 fields centered on Galactic 
cirrus clouds from which well matched on- and off- 
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cloud samples arc constructed. The main points of 
the paper are summarized below: 

• A CCD-based variant of the 75-yr old Wolf di- 
agram method is used to measure the effective 
extinction in different bands for the four cir- 
rus clouds. Star counts-based UBVRI extinc- 
tion curves are presented for two clouds, RCrA 
and PV 1. The steepness of the PV 1 extinction 
curve implies Ry = Ay / E(B — V) <J 2, which is 
unusually low compared to the standard Galac- 
tic value of Ry = 3.1. 

• A comparison of normalized cumulative stellar 
color distributions on and off the cirrus cloud 
yields color excesses, E(m\ — V) for m\ = 
UBRI. The reddening laws for all four cirrus 
clouds are steep, indicative of Ry values which 
are significantly lower than the standard value. 

• The optical depth of the RCrA and PV 1 clouds 
is also determined from the shift of the stellar 
distribution in U—B vs B-V and B—VvsV—I 
color-color diagrams. This method depends on 
an assumed value of Rv which determines the 
length, and to a lesser extent, direction of the 
unit reddening vector. Sample incompleteness 
at the faint end is found to be unimportant for 
the color-color analysis. 

• The four cirrus clouds in this study appear to 
favor unsually low values of Rv (<^ 2), com- 
parable to the lowest values found in previous 
studies. This fits in with the clouds' low optical 
depth (Ay < 1) and the general trend of in- 
creasing Ry with increasing optical depth (and 
vice versa) that is observed for the interstellar 
medium of the Galaxy. The low Ry values in 
optically thin cirrus clouds is suggestive of a rel- 
ative overabundance of small dust grains. 

• Constraints on cloud distances are obtained by 
studying the cumulative star counts versus ap- 
parent magnitude relations on and off the cirrus 
in the context of the Bahcall-Soneira star count 
model of the Galaxy. The RCrA cloud is esti- 
mated to lie at a distance of 550-750 pc, while 
the distance to the PV 1 cloud is < 1 kpc. 

• The effects of a non-uniform dust column den- 
sity are illustrated through a simple two-phase 
model: a low extinction sheet with embedded 



high density clumps. A clumpy medium tends 
to flatten (i.e., decrease the degree of wave- 
length dependence of) the star counts-based 
effective extinction curve. It also causes the 
counts-based Ay estimate to be higher than es- 
timates based on color excess and color-color 
data and lower than the estimate based on far 
infrared flux. The RCrA Wolf diagram and stel- 
lar color data are readily explained in terms of a 
non-uniform medium in which 20% of the area 
is covered by dense clumps. 

• Estimates of the dust optical depth derived 
from UBVRI extinction and reddening mea- 
surements are in good agreement with those 
based on the 100 /im flux. The values of the 
Ay /-Fioo/mii ratios for the four clouds are con- 
sistent with the ratios found in other studies of 
diffuse Galactic dust clouds and follow the ex- 
pected trends in dust temperature. 

• It should be possible to apply the photometric 
techniques described in this paper to the vast 
Sloan Digital Sky Survey data set in order to 
carry out an extensive study of dust properties 
in the diffuse interstellar medium and to con- 
struct an extinction map. 
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Fig. 1. — Greyscale (negative) representation of the 
V-band CCD images of the four cirrus cloud fields. 
The bold white contours indicate the on-cloud ("onl" 
and "on2") and off-cloud ("off") regions, demarcated 
on the basis of 100 /zm brightness. North is up and 
east is to the left. The scale of each image is indicated; 
the PV 1 image covers 1.8 deg 2 , while each of the other 
three images covers about 1 deg 2 . The optical surface 
brightness enhancement visible in the on-cloud region 
(relative to the off-cloud region) is due to reprocessing 
(mostly scattering) of ambient starlight by interstellar 
dust grains. The RCrA field is at the lowest Galactic 
latitude of the four fields (\b\ = 17°) and consequently 
contains the highest surface density of Galactic field 
stars. Areas around bright, badly saturated stars are 
excluded from data analysis. 

Fig. 2. — Cumulative star counts versus calibrated 
apparent L/-band magnitude (Wolf diagrams) in the 
four cirrus cloud fields. The y-axis numbers refer to 
the observed counts in the "off" region (solid black 
line); the "onl" (solid green line) and "on2" (solid red 
line) counts have been normalized to match the area 
of the "off" region (the "onl" counts are not shown 
for Paleyl and Paley3). The total number of satu- 
rated stars in each field and in each band (for which 
photometry is unavailable) is added to the bright 
end of the cumulative distribution. As expected, the 
"on2" (region with highest 100 /xm brightness) shifts 
are larger than the "onl" shifts. The dotted blue lines 
indicate Bahcall-Soneira star count models for the 
"off" and "on2" re gion s (upper and lower curve, re- 
spectively; see Sec. |J and Fig. |) of RCrA and PV 1, 
and for the "off" region of Paleyl and Paley3. The 
overall shape of the observed cumulative star counts 
(roughly a power law, with slight curvature) agrees 
well with that of the model relations, except for the 
flattening due to sample incompleteness at the faint 
end. The apparent U magnitude at which the sam- 
ple is 50% complete, U50, is indicated by a vertical 
dashed line for each field. 

Fig. 3. — Same as Figure || for the RCrA field in the 
BVRI bands, with Bahcall-Soneira star count models 
shown only for the B and V bands. Note the general 
trend of decreasing shifts between on-cloud and off- 
cloud star count relations as one moves towards longer 
wavelengths (m\ = U—*I). 

Fig. 4. — The apparent magnitude shift between 011- 
and off-cloud cumulative star counts, J 4= ounts ; plotted 



as a function of calibrated apparent magnitude m\ 
in the UBVRI bands (light blue, dark blue, green, 
red, and black lines, respectively) for the four cirrus 
clouds. The "onl" regions of Paleyl and Paley3 are 
not shown as the shifts are only marginally signifi- 
cant. As expected, the shifts decrease systematically 
with increasing wavelength from U to I for a given 
cloud/region and tend to be higher in "on2" than 
"onl". 

Fig. 5. — Estimate of the uncertainty in the founts 
values derived from star counts in the U and / bands 
(solid and dashed lines, respectively) for RCrA (field 
with highest star density) and Paley3. The ±ler er- 
ror band, plotted as a function of calibrated apparent 
magnitude m\ (U or /), is derived from the statisti- 
cal uncertainty in the star counts, projected onto the 
magnitude axis using the slope S\ of the log[V(< m)]- 
vs-m relation, combined with the photometric uncer- 
tainty. Poisson error dominates through most of the 
plotted apparent magnitude range; photometric er- 
ror and incompleteness cause the flaring of the error 
bands at the faint end. 

Fig. 6. — Effective extinction curve for RCrA and 
P V 1 based on ^4™ unts ; the median apparent mag- 
nitude (horizontal) offset between cumulative star 
counts in "off" versus "onl"/"on2" regions (squares 
and triangles, respectively, with la error bars), plot- 
ted as a function of inverse wavelength. The founts 
measurements for the other two clouds, Paleyl and 
Paley3, are not accurate enough for the computation 
of an extinction curve. The solid line shows the fidu- 
cial extinction curve from Cardelli et al. (1989) scaled 
to the weighted average value of the visual optical 
depth A^ nnts (see Table |) for the average Galactic 
value of R V =A V /E(B - V) = 3.1. For the RCrA cir- 
rus cloud, the fit is improved (relative to the uniform 
dust slab case, / = 1) if we assume that the inter- 
stellar medium is clumpy with (1 — /) = 0.1 of the 
cloud area covered by optically thick dense clumps 
of dust (dotted line); an identical fit can be made 
assuming Ry = 1.7 and (1 — /) = 0.2. For the 
PV 1 "on2" region, the fit is improved by adopting 
Rv — 1.7 (dashed line), which corresponds to the 
blue/ultraviolet portion of the dust extinction curve 
rising more steeply than for the standard extinction 
curve (Ry =3.1). 

Fig. 7. — Enlarged section of the t/-band Wolf dia- 
gram for the RCrA cloud (upper panel) and several 
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realizations of the Bahcall-Soneira star count model 
(center and lower panels). The star counts in the 
"off" , "on2" , and "onl" regions are indicated by bold 
solid, dotted, and thin solid lines, respectively, for 
the observations as well as the models. The U ~ 11- 
15 mag range shown represents the transition region 
between bright apparent magnitudes, for which the 
bulk of the stars are in front of the cirrus, and faint 
apparent magnitudes, for which most of the stars are 
in the background. Star count models are computed 
with and without an intervening slab of absorbing ma- 
terial, at distances of 0.4 kpc and 0.8 kpc that span 
the range of possibilities (center and lower panels, re- 
spectively), and with extinction corresponding to the 
measured A c jf unts = 0.89 mag in "on2" (left panel) 
and A c jf unts = 0.34 mag in "onl" (right panel). The 
la error bars in the center left panel are derived from 
an ensemble of model realizations each using the num- 
ber of stars appropriate for the actual RCrA "on2" 
area. 

Fig. 8. — Sample cumulative U — V (solid) and I —V 
(dashed) color distributions of stars in the "on2" 
(thin) and "off" (bold) regions of the RCrA and Pa- 
ley 3 fields. The horizontal bar indicates the median 
color excess shift measured in each case. An arbi- 
trary offset (in magnitudes) has been applied to each 
m\ — V (m\ — U or I) color scale in each field for 
convenience in plotting. 

Fig. 9. — Color excess, E(m\ — V) = A\ - A v 
(A = UBRI), derived from the median shift in the 
cumulative stellar color distribution (Fig. ||) of "onl" 
(squares) and "on2" (triangles) with respect to the 
"off" region, plotted as a function of inverse wave- 
length for the four cirrus fields. A standard Ry = 3.1 
Cardelli et al. (1989) reddening law scaled to the 
weighted average value of the visual optical depth 
^excess ( gee Table |j is plotted for the « onl » (dashed 

line) and "on2" region (solid line). The dotted line 
shows the best- fit Ry = 1.7 curve for the "on2" re- 
gion. Only the "on2" errorbars are shown for the sake 
of clarity. 

Fig. 10. — Color-color diagrams, U — B versus B — V 
and V — I versus B — V, for stars in the "off" (black), 
"onl" (green), and "on2" (red) regions of the RCrA 
and P V 1 fields (overplotted in this order) . A shift be- 
tween the on- and the off-cloud stellar distributions is 
visible, most clearly in the case of the RCrA UBV 
plot. Theoretical color-color loci from the Bertclli 



et al. (1994) stellar evolution models are shown as 
solid blue lines, averaged over all available ages for 
three different metallicities: Z = 0.0004, Z = 0.008, 
and Z = 0.05 (top to bottom). Reddening vectors 
corresponding to Ay = 2 mag are plotted for the 
standard Ry = 3.1 extinction law (red arrow) and 
for Ry = 1.7 (black arrow). A given offset in the 
color-color plane translates to a larger Ay for a stan- 
dard interstellar reddening law than for an Ry = 1.7 
law. 

Fig. 11. — Illustration of the shift of the on-cloud 
["onl" (dashed) and "on2" (solid)] stellar UBV color- 
color distribution with respect to the off-cloud distri- 
bution, measured along the Ry = 3.1 reddening vec- 
tor, for the RCrA field. The stellar distributions are 
compared in a series of 0.1 mag- wide slices parallel to 
the reddening vector, and A is the offset of each slice 
in the direction perpendicular to the reddening vec- 
tor relative to an arbitrary origin. The median offset 
between "on2"/"onl" and "off" stellar distributions 
along each slice is plotted, along with la error bars 
based on the width of the distribution and the number 
of stars. The dot-dashed horizontal line indicates the 
best- fit value of the visual optical depth, Ay BV , de- 
rived from the weighted mean of the offsets measured 
in the individual slices. 

Fig. 12. — The effect of dumpiness on the shape of the 
counts-based extinction curve, plotted in the normal- 
ized form: ^counts ^counts versus inverse wavelength. 
The calculation is based on a simple two-phase model 
for the cirrus cloud (see Sec. |6.2| ) comprised of a 
smooth sheet of absorbing material with optical depth 
Ay cct and dense, opaque clumps (A v ump — > oo). 
A one-parameter family of model counts-based ex- 
tinction curves is obtained by varying the ratio, 

^sheey^counte for ft fixed yame of ^counts = ^ tne 

overall offset in the 1^-band Wolf diagram. The curves 
shown are computed using the star count slopes S\ 
measured in the RCrA field and a standard Ry = 3.1 
extinction law. The quantity / is area covering factor 
of the low-optical depth, "sheet" portion of the cloud, 
and the listed values are based on ^4™ unts = l. The 
higher the degree of non-uniformity of the cloud, the 
flatter the shape of the counts-based effective extinc- 
tion curve; the limiting shape (for Ay cct /Ay unts — 
0) is determined only by the relative values of the 
UBVRI star count slopes Siambda. 
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Fig. 13. — Comparison of Ay estimates derived 
using various methods (left to right): scaled from 
100 /im brightness and far infrared color tempera- 
ture measurements (Sec. ||), cumulative star counts 
in the UBVRI bands (Sec. |^), cumulative color dis- 
tributions (Sec. 
diagrams (Sec. 



5.1), and UBV and BVI color-color 



5.2). Most error bars do not 



ceed the symbol sizes, and, for the sake of clarity, 
have not been plotted. There is good agreement 
amongst the different Ay estimates. The dotted line 
in panel 1 (RCrA "on2" ) shows Ay values for a two- 
phase clumpy dust model (Sec. 6/2; Fig. |lj) in which 
a fraction (1 — /) = 0.2 of the area is covered by 
optically-thick clumps; the range of ^4^P°^ m values 
(IR) indicated corresponds to Ay — 2.4 - 3.5 mag 
in these dense clumps. Open symbols indicate Ay 
estimates based on the standard Ry = 3.1 interstel- 
lar extinction law; filled symbols are based on a non- 
standard extinction law (Ry = 1.7) which appears 
to be favored by the color excess data in all fields. 
For Paleyl and Palcy3, only the "on2" estimates are 
reliable and these are plotted. 
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